Electron Phonon Superconductivity in LaNiOP 
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We report first principles calculations of the electronic structure, phonon dispersions and elec- 
tron phonon coupling of LaNiPO. These calculations show that this material can be explained as 
a conventional electron phonon superconductor in contrast to the FeAs based high temperature 
superconductors. 

PACS numbers: 74.25. Jb,74.25.Kc,74.70.Dd,71.18.+y 



The discovery of high temperature superconductiv- 
ity in a series of FeAs based compounds, prototype 
LaFeAs(0,F) by Kamihara and co-workersi has led to 
intense investigation of these phases and the discovery 
of a number of new phases with critical temperatures in 
some cases near and above In addition, recently, 

high temperature superconductivity with similar proper- 
ties has been found in simpler compounds based on FeAs 
sheets such as ThCr2Si2 structure type (Ba,K)Fe2As2-— & 
These Fe-based materials are commonly viewed as uncon- 
ventional (non-electron-phonon) superconductors based 
on the high values of T c , the proximity to magnetism 
and the fact that calculations of the electron phonon 
coupling show that it is far too weak to account for the 
superconductivity^ It should be emphasized that these 
are rich structure types and that there is no doubt much 
compositional space that remains to be explored!*^ 

In this regard, besides the FeAs based phases, su- 
perconductivity is also observed in Ni based materi- 
als, including LaNiPO ) 11 ' 12 pure, fluorine or Sr doped 
LaNiAsO^iii^ and LaNiBiO^ Interestingly, both Ni 
and Fc are ambient temperature elemental ferromagncts, 
and many Fe and Ni compounds show magnetism. In 
addition, it is interesting to note a possible connection to 
the cuprates. Specifically, the Ni compounds are based 
on square lattices of nominally Ni 2+ and superconduct 
with various dopings, including electron doping to yield 
a nominal d occupancy 8+5, while the high T c supercon- 
ductors are based on square lattices of Cu 2+ with hole 
doping for nominal d occupancy 9-<5, where 5 is the dop- 
ing level. 

Returning to the connection with the Fe-based com- 
pounds, both electronic structure calculation s] 17 ' 18 ' 19 ' 20 
and experiment show important differences between the 
Ni-based and Fe-based superconductors. In particular, 
the Fe-based materials show low carrier density with rel- 
atively small Fermi surfaces, high density of states and 
proximity to magnetism ; 8 ! 17 ' 21 ' 22 while the Ni compounds 
show large Fermi surfaces, lower density of states and are 
apparently further from magnetic instabilities. This sug- 
gests that the Ni-based compounds may be a different 
class of superconductors, perhaps, considering the lower 
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FIG. 1: Calculated LDA band structure of LaNiPO using the 
experimental lattice parameter and relaxed internal coordi- 
nates, ZLa=0.1506 and zp =0.6210. 



T c observed so far in this group, conventional electron- 
phonon materials. On the other hand, it may be noted 
that the kinds of measurements and calculations done 
so far would have given qualitatively similar results for 
cuprate superconductors in the the optimal and over- 
doped regions - i.e. moderate density of states, high 
carrier density metals, apparently far from magnetism. 

Here we report details of our previous electronic struc- 
ture calculations as well as calculations of the phonon 
dispersions and electron-phonon coupling for LaNiPO. 
The calculations show that unlike the Fe-based materials, 
superconductivity in this Ni-based compound is readily 
explained by the standard electron-phonon mechanism. 
This means that the superconductivity of LaNiPO is not 
related to that of the FeAs based materials. 

Our electronic structure was discussed briefly 
previously^ 7 - and is similar to that reported by Zhang 
and co-workers for the same compound^ It also shows 
similarities to the electronic structures found for the 
other Ni-based compounds . 19 ' 20 Our electronic structure 
calculations were performed within the local density 
approximation (LDA) using the general potential lin- 



FIG. 2: (color online) Calculated LDA density of states of La- 
NiPO on a per formula unit both spins basis. The projections 
are onto the LAPW spheres. 





FIG. 3: (color online) Calculated LDA density of states of 
LaFePO on a per formula unit both spins basis. The projec- 
tions are onto the LAPW spheres. The Fe sphere radius was 
2.1 ao, which is equal to the Ni sphere radius used in Fig. [5] 



earized augmented planewave (LAPW) method 2 ^ as for 
LaFeAsO (Ref. LAPW sphere radii of 2.2 o , 2.1 a , 
2.1 ao and 1.6 a were used for La, Ni, P, and O, respec- 
tively. In particular we used the experimentally reported 
tetragonal lattice parameters (a=4.046lA, c=8.100A)ii 
and relaxed the internal coordinates which correspond 
to the La and P heights. We obtain ZL a =0.1506 and 
zp=0.6210, which are close to the reported experimental 
values of z La =0.1531 and z P =0.6260,Ii and z La =0.1519 
and zp=0.6257»i£ This is different from the Fe-based 
compounds where pnictogen heights significantly lower 
than the reported experimental values arc generally 
obtained^i^ 

The calculated band structure and electronic density 
of states (DOS) are shown in Figs. Q] and [21 respec- 
tively. The DOS for LaFePO, calculated in the same 
way is shown in Fig. [3] for comparison. The Fermi sur- 
face of LaNiPO is given in Fig. 0) As may be seen, these 
are very different in LaNiPO from those in LaFeAsO. 




FIG. 4: (color online) Calculated LDA Fermi surface of La- 
NiPO shaded by velocity. The top panel shows a view along 
the tetragonal axis while the bottom panel shows a tilted view. 



This is as might be expected from the different electron 
count. LaNiPO and LaFeAsO have a similar structure 
to the d bands, and in particular LaNiPO should be de- 
scribed as Ni 2+ ions on a square lattice with direct hop- 
ping as well as a P induced crystal field. In fact, the 
projections of the density of states shows that there is 
a greater degree of covalcncy between Ni and P, than 
between Fe and As or Fe and P in the corresponding 
Fe-based materials. There remains a pseudogap at a d 
electron count of 6, however Ni 2+ has 8 electrons. This 
places Ep well inside the upper manifold of d states. 
In this energy range the bands are derived from Ni d 
states hybridized with P p states. The bands at Ep are 
more dispersive than in the Fe compounds where Ep is 
lower. This leads to lower density of states with higher 
in plane Fermi velocity, N(Ep)=lAl cV _1 per formula 
unit, v xx =:3.75xl0 7 cm/s, v Z2 =0.39xl0 7 cm/s. The lower 
N(Ep) puts the Ni based compound further from mag- 
netism than the Fe-based materials, as was discussed*^ 
Furthermore, this compound is quite two dimensional in 
the sense that there is no 3D Fermi surface sheet. Based 
on the anisotropy of the Fermi velocity, the ratio of in- 
plane to c-axis conductivity for the Ni compound is ~ 100 
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FIG. 5: Calculated phonon dispersions of LaNiPO as obtained 
within linear response. 




FIG. 6: Calculated phonon density of states and electron 
phonon spectral functional a 2 F(to). 



assuming isotropic scattering. The Fermi surface may de- 
scribed as consisting of two large ellipsoidal cross-section 
electron cylinders around M, and a large hole section 
around T. This hole section intersects the zone boundary 
near X leading to an electron section around X. 

The phonon and electron phonon calculations were 
performed in linear response, with the quantum espresso 
code2£ and ultrasoft pseudopotentials and the general- 
ized gradient approximation of Perdcw. Burke and Ernz- 
erhof (PBE),— similar to our previous calculations for 
LaFcAsO.— We did convergence tests for the basis set size, 
the planewave expansion of the charge density, the tem- 
perature broadening, and the Brillouin zone sampling. 
An 8x8x4 grid was used for the zone integration in the 
phonon calculations, while a more dense 32x32x8 grid 
was used for the zone integration in the electron-phonon 
coupling calculation. The basis set cut-off for the wave- 
functions was 40 Ry, while a 400 Ry cut-off was used for 



the charge density. 

The calculated phonon dispersions of LaNiPO are 
shown in Fig. [5] The corresponding phonon density 
of states and Eliashberg spectral function a 2 F(cu) are 
shown in Fig. [5] The phonon dispersions show a set of 
12 phonon bands extending up to ~ 200 cm -1 , separated 
by a gap from 12 higher frequency bands extending up 
to ~ 500 cm -1 (note that there are 24 phonon branches 
since there are two formula units per primitive cell). The 
higher frequency manifold is derived mainly from O and 
P motions. Within this upper manifold the P contribu- 
tion is mainly below ~ 300 cm" 1 , while the dispersive 
modes above 300 cm -1 are mainly O derived. The lower 
manifold from to 200 cm -1 consists of the acoustic 
modes and modes of mixed, but mainly metal character. 

We obtain a value of the electron phonon coupling 
A=0.58 with logarithmically averaged frequency 0^=113 
cm" 1 . Relative to the phonon density of states, the spec- 
tral function is enhanced for the lower frequency metal 
modes, which have strong in-plane Ni character, and for 
the modes at the bottom of the upper manifold, which 
have strong P character. Thus in spite of the lower 
electronic density of states, we obtain stronger electron 
phonon coupling as compared with LaFeAsO, where A ~ 
0.2ii& Inserting these numbers into the simplified Allen- 
Dynes formula, 



k R T r 



in 



1.2 



cxp 



1.04(1 + A) 
A-^*(1 + 0.62A) 



(1) 



with ordinary values of the Coulomb parameter fi* 
yields values in reasonable accord with experiment. 
Specifically, for fi* =0.12 we obtain T C =2.6K, which is 



3K, 



or 



in accord with the experimental value T c 
T C =A.2KH 

In conclusion we find that LaNiPO has a conventional 
superconducting state which arises from band metal with 
moderate density of states and intermediate electron 
phonon coupling. This is in contrast to e.g. LaFeAs(0,F) 
which is a high density of states material, near mag- 
netism, and with weak electron phonon coupling that 
can in no way explain the superconductivity. This leav- 
ing aside structure and chemistry, the superconductivity 
of LaNiPO and presumably the rest of the Ni-based oxyp- 
nictides is unrelated to that of the Fe-based materials. 

We are grateful for discussions with I.I. Mazin, D.G. 
Mandrus, B.C. Sales, R. Jin, and M. Fornari and support 
from DOE, Division of Materials Sciences and Engineer- 
ing. 
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